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We propose to use a quantum ratchet to transport quantum information in a chain of atoms 
trapped in an optical superlattice. The quantum ratchet is created by a continuous modulation of 
the optical superlattice which is periodic in time and in space. Though there is zero average force 
acting on the atoms, we show that indeed the ratchet effect permits atoms on even and odd sites 
to move along opposite directions. By loading the optical lattice with two-level bosonic atoms, this 
scheme permits to perfectly transport a qubit or entangled state imprinted in one or more atoms to 
any desired position in the lattice. From the quantum computation point of view, the transport is 
achieved by a smooth concatenation of perfect swap gates. We analyze setups with noninteracting 
and interacting particles and in the latter case we use the tools of optimal control to design optimal 
modulations. We also discuss the feasibility of this method in current experiments. 



I. INTRODUCTION 

By a ratchet effect one usually refers to the existence 
of directed transport in a system in which there is no 
net bias force. Ratchets have been traditionally found 
in the study of dissipative systems where external 
fluctuations causing Brownian motion cooperate with a 
periodic force to bias transport, a mechanism that is in 
the basis of some biological motors [2[ . Ratchets can also 
appear without dissipation. These Hamiltonian or con- 
servative ratchets are interesting as they can be extended 
to quantum mechanical systems. In this context one finds 
studies that relate the existence of transport to classical 
properties of the model, such as some asymmetries of the 
external force J3| or mixing of chaotic and regular phase 
space regions [J]. 

Ultracold atoms offer an ideal arena to test this phe- 
nomenology, as can be seen from the experiments im- 
plementating both d issip ative and conservative ratchets 
E 1 0, E S M EL El 0, El El • These experiments 
rely on the force imparted by near or far from resonance 
laser beams which act on the atoms over short periods 
of time. These flashing potentials implement variants of 
the ^-kicked rotor model and lead to phenomena such as 
dynamical localization. On the theoretical side we must 
remark the discovery of directed transport on quantum 
mechanical systems whose classical counterpart is com- 
pletely chaotic [H, E3, El, Ei H3 . 

The present work aims at exploiting the fact that quan- 
tum ratchets can be used to transport quantum states 
and thus distribute entanglement between distant nodes, 
a basic task for quantum information and computation 
purposes [2l|. Within this context, we find two dis- 
tinct research directions. On the one hand, entanglement 
swapping and quantum teleportati qn |22l c an be used to 
build efficient quantum repeaters [23l l24j for long dis- 
tance communication. On the other hand, on a smaller 
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FIG. 1: (Color online)An optical superlattice arises from a 
combination of two potentials with different periods. By mod- 
ulating the depths and displacements of these potentials we 
can raise and lower the tunneling rates between odd and even 
pairs of sites (red arrows). This way, by means of perfect 
swaps, the state of a particle can be transported along the 
lattice. 



scale, a relevant effort has been recently devoted to study 
quantum state transfer using static local Hamiltonians 
which act on chains of spins, harmonic oscillators, or 
bosons in optical lattices J2p, IM III III IM S3, IM SI, 
HI, M, SE SI |U Compared to these works, 
time-dependent (non-adiabatic) models, such as the one 
presented here, provide more degrees of freedom to opti- 
mize the efficiency, speed and robustness of the quantum 
transport. 

Moreover, the scheme we propose here for quantum 
state transport can be seen as a quantum ratchet in- 
duced without breaking translational and time reversal 
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symmetries in the forces. The transport is just triggered 
by the symmetry breaking of the initial state. We illus- 
trate this with a particular and simple implementation of 
quantum state transport based on two-level bosonic ul- 
tracold atoms in optical superlattices. More specifically 
the transport is achieved as follows. All the atoms are 
initially prepared in the, say, down state. At a given 
site (write port) an arbitrary qubit state is imprinted 
in the atom. The trapping potential is then modulated 
smoothly and periodically in time and in space. Depend- 
ing on the initial asymmetric position of the atom, being 
in an even or odd lattice site, the qubit state is trans- 
ported rightwards or leftwards to any desired lattice site 
(read port) [Fig. Q] without precise individual addressing 
[53j | . We consider both free and interacting particles. To 
our knowledge, this is the first time many-body interac- 
tions are taken into account in Hamiltonian ratchets. 

From the implementation point of view, this work is 
inspired by recent advances in the control and manipula- 
tion of optical superlattices [HI, H2, EE E3, El ■ Present 
experiments can create periodic potentials such as the 
ones depicted in Fig. [TJ achieving a great control on the 
time-dependence of the potential heights, and even being 
able to measure the number of atoms on the even and odd 
sites [H, El| ■ These tools suffice for the protocol devised 
in this work. 

The outline of the paper is as follows. In Sect. HT1 we 
introduce our system, made of a chain of atoms in an op- 
tical superlattice. We state the goal of this work, which is 
to transfer a quantum state between two arbitrary lattice 
sites using a translationally invariant modulation of the 
superlattice, and we present the Bosc-Hubbard Hamilto- 
nian that models the dynamics of the atoms. In Sect. IIIII 
we study the case of noninteracting particles and find 
that it is possible to induce directed transport with arbi- 
trary speed. In Sect. [IV] we discuss a more realistic setup 
in which atoms interact. We design a ratchet Hamil- 
tonian with the help of optimal quantum control, and 
find that we can still induce perfect quantum transport 
with a speed limited by the interaction strength. The 
details of the optimal quantum control technique are left 
to Sect.|VTl In Sect. [V] we study the optical superlattice 
modulation in more detail, discussing the experimental 
challenges for implementing these quantum ratchets in 
current experiments. Finally, in Sect. IViTl we summarize 
our results and comment on possible extensions. 



II. OPTICAL SUPERLATTICES 



or maxima of intensity. The optical superlattice arises 
from a combination of two potentials with different peri- 
ods. In the simplest case the superlattice potential can 
be written as follows 



V{x,y,z,t) = V x (t) cos 2 (kx) + V 2 (t) cos 2 (2kx 
+ V± [cos 2 (ky) + cos 2 (kz)). 



(1) 



Here V x (t) and V± are the strengths of the lattice along 
the main axis and transversely to it, and k — 2ir/\ is the 
momentum of the photons. In order to have a configu- 
ration of decoupled ID lattices [47], the transverse po- 
tential V± has to be much larger than the recoil energy 
E r = ti 2 k 2 /2m to prevent tunneling between ID lattices. 
On top of this potential and along the main axis we find 
another lattice of strength V^t), with half the period of 
the original one and a possible dephasing <f> £ Z x ir/2. 
This yields to the optical superlattice that creates the 
double-well structure shown in Fig. [T] 

We are interested in the regime of deep lattices in 
which atoms are confined to the lowest Bloch band of the 
periodic potential. We assume the atoms to be bosons 
with two internal states. Under those circumstances the 
dynamics of the atoms can be modeled using a Bose- 
Hubbard Hamiltonian [46j which, for a single ID tube of 
lattices, has the form 
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(2) 



Here, a a i and a} ai are Fock operators that annihilate and 
create particles on the z-th site of the lattice and in one 
of two internal states, a E {t, !}■ The parameter U is the 
on-site interaction between atoms, which we assume to 
be spin independent, and J t is the tunneling amplitude, 
which may vary from site to site. 

In order to describe the dynamics of the atoms, Eq. (|2|) 
has to satisfy two conditions. First of all, the tunnel- 
ing Jj between different wells has to be small, so that 
there effectively exists a low energy band formed by the 
localized states on each site. This implies that the su- 
perlattice cannot drop the barrier between pairs of sites 
too much [Fig. [Tj, for otherwise we would have to use 
a multiband model. The second condition is that the 
on-site interaction U must be small compared to the 
energy gap between Bloch bands. Both condition are 
easily met in current experiments with optical superlat- 
tices [il|, El, El, El| , and have been taken into account 
throughout this work [See Sect. W\. 



A. The model 

In this work we consider the setting of ultracold neu- 
tral atoms trapped in an optical superlattice. An optical 
lattice is a standing wave of coherent off-resonance light 
created by the interference of two or more laser beams. 
As explained in Ref . [46] , such standing wave behaves as a 
periodic potential that confines the atom in the minima 



B. The goal: quantum communication 

Our objective is as follows. We initially prepare all the 
two-level atoms in the same state, say the down state 
||). We imprint an arbitrary quantum state on a partic- 
ular site of the lattice, the write port, in such a way that 
the atom ends up in a superposition of both spin states, 
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a IT) +/5 |4)- The goal is to create, by modulating the en- 
ergy barriers of the optical superlattice, a ratchet effect 
that will transport this state a predetermined distance, 
to another site, the read port in Fig. [TJ The modula- 
tion of the superlattice will not break the translational 
invariance and it will also be periodic in time (with a 
period equal to 2T) , alternating the roles of the hopping 
between even (J2„ = J 2) and odd (J2n+i = J\) sites 

Ji, 2 (t + T) = J 2 ,i(*). (3) 

The Hamiltonian in Eq. (J2]), being translationaly in- 
variant in time and in space, does not exert any force on 
the atoms, as expected from a quantum ratchet. This 
can be verified by computing the integral of the forces 
acting on an atom in a well of the lattice. In particular, 
the following average is zero 

r 2T r^k a] 

/ dt dx—V(x,y,z,t) = 0. (4) 
Jo Jo dx 

Note however that our superlattice does not have a per- 
manent asymmetry on the unit cell [0, 7r/2fc), unlike other 
Hamiltonian flashed ratchets with saw-tooth profiles Q. 

While the average force acting on the atoms is zero, 
we will show that there is still transport. The reason is 
the asymmetry of the initial state 

\m) - («a{ in + (3a\ in ) [J aljvac), (5) 

and the fact that our time-dependent Hamiltonian will 
perform perfect swaps between neighboring lattice sites. 
All of this leads, after n swaps of duration T, to a prop- 
agation of the imprinted qubit rightwards or leftwards 

|V(nT)) = (aa] out + (3a\ J J] aljvac), (6) 

depending on the starting site. Note also that by design, 
we will achieve the transport of the quantum state with- 
out actually moving any atom, thus leaving the number 
of atoms per site invariant. 

Finally, let us mention that this scheme can also be 
used to distribute entanglement in the same fashion as 
in (27j . Namely by imprinting a maximally entangled 
state in two neighboring lattice sites. Since atoms in odd 
and even sites move on opposite directions, the atoms 
will depart from each other and entanglement will get 
distributed between arbitrarily distant lattice sites. 

III. NONINTERACTING CASE 

We start with the case of nonintcracting particles 
U = 0. This case is relatively easy to study as all parti- 
cles can be treated independently and therefore we can 
focus on the dynamics of a single particle that starts from 
different sites, a problem which is integrable. The basis 



of states is now denoted by \j) :— aj|vac), where we no 
longer care about the internal state of the particle. Our 
only concern is now to ensure that the state with a par- 
ticle on the write port \j = in) is mapped after a given 
time, t r , to a state with that particle on site \j = out), 
the reading port, in a deterministic fashion. Note also 
that since particles do not interact we do not have to 
consider the dynamics of particles on other lattice sites. 
Nevertheless, as we show below, it is possible to look for 
solutions that do not give energy to these other particles 
and leave the density profile invariant. If the number of 
atoms per lattice site is constant, we can then talk about 
transport of the quantum state and not of particles them- 
selves [Fig. [T]. 

As described before, we will adopt two restrictions. 
The first one is that the couplings do not break the trans- 
lational invariance but are modulated in time, with the 
roles of J2 and J\ being exchanged after a time T, as from 
Eq. ([3]). The second one is to assume that the tunneling 
can be brought to zero, so that 

Ji(t)>0, J 2 (t) = 0, te [0,T), 

JiAt + T) := J 2 ,i(t). (7) 

During the time in which Ji = we have (jx = 1) 

where a x is a Pauli matrix and 9{t) = J Q </i(t)g?t. The 
simplest solution of this kind is a square signal, J\ (i) = J, 
as shown in Fig. O). This means that for T = 7r/(2J) we 
achieve a perfect swap between even and odd sites up 
to an irrelevant global phase. In practice the tunnelings 
Ji and J2 will evolve smoothly and require some time to 
reach a nonzero value. In that case the time for perfect 
switch will be given by 0{T) — it/2. 

It is possible to construct arbitrarily fast solutions of 
also arbitrary smoothness. Note, however, that the value 
of the average hopping is inversely proportional to the 
modulation period, J cx 1/T, so that a fast solution may 
require a lattice with an unrcalistically large hopping. 

IV. INTERACTING CASE 

The case of interacting particles is more realistic but 
also more difficult. We can no longer regard the parti- 
cles as independent and collisions can affect the phase of 
the transported state, as defined by Eq. ([5]). We will 
still look for simple solutions that concatenate a swap 
dynamics on pairs of lattice sites, with the restrictions 
introduced in Sec. Ill Bl 

A. A two-level problem 

Let us repeat the calculation of Sect. QTI] but now tak- 
ing into account the interaction between particles. Since 
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FIG. 2: (Color online) (a) Model case of two lattice sites 
disconnected from the rest because J2 = 0. By fixing the 
hopping between wells to a precise value J for a time T, it is 
possible to swap the atoms. Note that the doubly occupied 
sites have more energy due to the interaction, U, and a possi- 
ble inhomogeneity of the lattice potential, A. (b) Combining 
the solution for a pair of sites we obtain one possible mod- 
ulation of the hoppings Ji (solid) and J2 (dashed), namely 
Ji and J2 are square waves in antiphase. This produces the 
quantum transport, (c) Values of J and T for which perfect 
transport is achieved. Each circle is a solution; the solid line 
joins the solutions with smallest hopping and the dashed line 
is the solution for noninteracting particles. 



J2 = we focus on the double-well problem with only 
two distinguishable particles, i.e. in two different inter- 
nal states. Instead of using second quantization we de- 
note their state as \i)\j), where i,j € {0,1} is the well 



in which the particle resides. The connections between 
states are depicted in Fig. [5^,. Using Pauli operators the 
effective Hamiltonian becomes y54f 



H = -J(a x ® I + I <g> a x ) + ^(a z 



1 



(9) 



We can further restrict our problem to the only states 
that participate on the dynamics 
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(|00) + |11». 



(10) 



Using this basis we obtain the effective Hamiltonian 


H(t) = ( 





(11) 



where \ip~) is shown to be a dark state and the other two 
are coupled by a two-level Hamiltonian 2J(t)a x + Ua z . 
Our problem is thus to find a hopping J(t) such that after 
a time T the states above have experienced the following 
transformation 

hn - hn, i^ + > - -i^ + >, i^ + > -> e^i0+>, (12) 

where the phase v is unimportant for our purposes. 

B. The square signal revisited 

We will first investigate solutions which are piecewise 
constant, with J\{t) — J, for t S [0, T). For a fixed 
hopping our Hamiltonian has two nonzero eigenvalues 



E± = 



C/±V16J 2 + L/ 2 



(13) 



that contribute to the evolution of the symmetric states. 
The perfect swap between the atoms takes place at a 
time T such that the symmetric state \ip + ) changes sign. 
As can be written as a superposition of the two 

symmetric eigenstates with energies E±, this condition 
is satisfied when E±T — ±(2n± + 1)71", where n+ and n_ 
are arbitrary integers. Defining x := {2n + + 1)/ (2?i_ + 1) 
we obtain 



J \/X 



(14) 



From this value we can compute the energies E± and 
the time T. As shown in Fig. [2t the minimal value of 
J is reached for either n + or n_ equal to and, unlike 
in the noninteracting case, there exists a minimal swap 
time given by T — 2n/U. Indeed, using the tools in [48| 
one can prove that the constant hopping is the fastest 
solution and that the gate cannot be performed faster 
than this time. 
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FIG. 3: (Color online) (a) Hopping and (b) fidelity of the gate 
(F = |(01|t/(i)|10)| 2 ) with the perfect transport for half a pe- 
riod [0,T) during which Ji(t) = J(t) and J 2 = 0. We find 
two solutions, one for T = 2n/U (solid) and another one for 
T = 4n/U (dashed), with three and two modes, respectively. 
In (c) we plot the corresponding modulations of the superlat- 
tice. 

C. Smooth solutions: optimal quantum control 

While optimal, the square signal that we have con- 
sidered before is probably unrealistic, as in experiments 
the hopping will smoothly increase and decrease as the 
tunneling barriers are changed. For that reason we have 
investigated other solutions with continuous derivatives 
using the tools of optimal quantum control. 

The details of the method are left for Sect. IVI1 but let 
us sketch the procedure. The first step is to parametrize 
the hopping as a linear combination of some functions 

J(t)=Y, C nfn(t)>0 (15) 
n 

where for simplicity we use trigonometric functions 
f n (t) = sin 2 (nnt/T) that increase and decrease smoothly 
to zero. Since the hopping is positive we have a first 
restriction, c„ > 0. The second restriction is that the 
fidelity of the swap procedure has to be one. Both con- 
straints are imposed to the problem of minimizing the 




FIG. 4: (Color online)On the upper figure we plot a solution 
of Ji(t) (solid) and J 2 (i) (dashed) for a perfect transport of 
the qubit. In the lower figure we plot the average position of 
the qubit on the lattice as transported by these modulations. 
These plots have been computed using T = 4ir/U. 

average strength of the hopping, given by E = J2 n l c «l • 
This problem is solved numerically with MATLAB's op- 
timization toolbox 49] using as aid the derivatives com- 
puted by means of perturbation theory [Sect. I VI] . 

In Fig. [3] we show two instances of the problem, one 
optimized for three modes (c„ = for n > 3) and a 
duration of T = 2n/U, the other one for twice the time 
and two modes. As it can be appreciated in the picture, 
we reach perfect fidelity in a rather smooth and robust 
manner, so that errors in the timing of the gate will not 
affect the process significantly. 

To further relate these solutions to the notion of a 
ratchet let us look at Fig. [U There we plot the full time 
evolution of the imprinted qubit state for 6 periods of an 
optimal modulation. The position may oscillate between 
pairs wells, but there is always a net average transport. 



D. Dealing with holes 

The two solutions derived above, that is the piecewise 
constant and the smooth ones, are designed to induce 
transport on a chain of particles. However, in practice 
such a chain will have some endpoints or particles that 
stand near an empty site. We can then have three sce- 
narios: (i) that a particle standing near a hole ends up in 
the original site at time t = T; (ii) that the particle and 
the hole are swapped and, more generally, (iii) that the 
particle and the hole end up in some coherent superpo- 
sition of being on each site. Out of these processes, only 
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FIG. 5: (Color online) Parameters of the double well potential 
as a function of the superlattice modulation, AV = (V x + 
V2)/2. We plot the effective hopping between wells, J (dash- 
dot, right axis), the energy gap to higher bands, AE (dashed), 
and the on-site interaction, u (solid) . Everything is expressed 
in units of the recoil energy, E r . 



the latter will affect the evolution of the state we want to 
transport, since there is a small probability that it gets 
reflected. 

Holes only have a disturbing effect on the transport if 
U 7^ 0, since in the noninteracting case surrounding par- 
ticles are equivalent to holes. Note however that we can 
effectively eliminate the scenario (iii) if we impose that 
the hopping on half a period leaves the particle invariant 



[ J{r)dT = 2tt x Z, 
Jo 



(16) 



While this restriction cannot always be achieved for the 
piecewise constant profile, it can be easily incorporated 
to our optimal control toolbox. 



EXPERIMENTAL IMPLEMENTATION 



and V2 [See Eq. {T])] . In order to relate these two quanti- 
ties we have have performed a band structure calculation. 
We have focused on the case J2 (t) ~ 0, which corresponds 
to the first half-period, t £ [0,T). This hopping can be 
suppressed by making (j) — and V x + V2 = 70E r , which 
is large enough to effectively suppress hopping every sec- 
ond site. Given this constraint, any modulation of the 
lattices depends on a single quantitiy, AV(t) > 0, such 
that V x = AV(t) and V 2 = 7QE r -AV(t). It thus remains 
to relate AV(t) and J\{t). This is done for each possi- 
ble value of AV, by computing numerically the ground 
state wavefunction of a particle in a double- well, w(x), 
and the first excited state. The energy difference be- 
tween these states is proportional to the effective tunnel- 
ing Ji(t), while the on-site interaction energy becomes 



U_ 



(17) 



where a s is the scattering length between atoms, a the 
period of the superlattice and u = J \w(x)\ 4 . 

The results are shown in Fig. [5] Since a s /a is very 
small, the on-site interaction energy is smaller than the 
energy gap between Bloch bands and the Bose-Hubbard 
model is therefore valid throughout the evolution (|2|). 
Furthermore, u does not change very much as a func- 
tion of the modulation AV, while J decays exponentially 
fast with AV. After an appropiate fit of these quanti- 
ties, one may convert the time dependence of J\(t) from 
Fig. [3^ into the associated modulations of the superlat- 
tice AV(t), which are shown in Fig.[3j:. 



VI. COHERENT CONTROL 

In this section we present the tools that we have used 
to optimize the modulation of the hopping, as they differ 
significantly from what is the standard approach in opti- 
mal quantum control based on a Lagrangian formulation 



All the protocols that we have designed can be imple- 
mented in current experiments with optical superlattices 
[U S3, |H, Q HH. The implementation should begin 
by loading the superlattice with approximately one atom 
per site, all of them in the same internal state. The next 
step is to rotate the state of some atoms using either 
magnetic field gradients, coherent light or a clever com- 
bination of both [50| , in order to imprint the quantum 
state we wish to transport. By modulating in time the 
intensity of the laser beams that participate in the opti- 
cal lattice one can then achieve a modulation of J\ (t) and 
J>2,{t) that corresponds to the solutions studied above. Af- 
ter an appropiate time one may retrieve the qubit state 
by measuring the lattice sites on which it is expected to 
arrive. 

In an experiment one does not directly control the hop- 
pings Ji(t) and hit), but rather the lattice strengths V x 



A. Objective function 

Let us formulate our problem: we have a Hamilto- 
nian, Hit; x) that depends both on time and on some 
additional parameters, x\ . . . xm- Our goal is to find an 
optimal set of parameters of the Hamiltonian such that 
the evolution of a number of states is as equivalent as 
possible to a transformation given by a specific unitary 
U g . 

Stated in a more concrete way, the evolution of an ar- 
bitrary state IV'(O)) is given by \ip(t)) = x)|^(0)), 
where the unitary is a solution of the Schrodinger equa- 
tion 

i^-Uit;x) =H{t\x)U{t;x) (18) 
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with initial condition U(0;x) = I [55|]. Our goal is to 
maximize the fidelity (defined below) of any evolved state 
U (t; x)\ip n ) with the desired transformed state U g \tj} n ). 
There are many ways to measure the accuracy of the 
transformation. A strict and simple objective function is 
the fidelity 

1 1 d 

F = -Re{tr[[/t(7(T;x)]} = -Re Y,{^Wp{T; x)\j>n), 

n=X 

(19) 

where the {IVVi)} form an orthonormal basis of the d- 
dimensional Hilbert space H, where we want to control 
the evolution (56[. This function is bounded by F £ 
[— 1, 1] and achieves the maximum value for the perfect 
transformation, 



F = 1 U(T;x) 



U g . 



(20) 



The question is thus, how do we maximize Fl A natu- 
ral way is to compute (when possible) the derivative of 
F with respect to the parameters x, i.e. dF/dxi, since 
the gradient itself provides a direction along which the 
fidelity is increased. Indeed, given this derivative there 
are multiple optimization algorithms that would allow us 
to compute the optimal control. 



B. Formal gradient 

In order to obtain the gradient of F with respect to 
the parameters x, we straightforwardly obtain from (Q2 



Ipn 



(21) 



and same initial condition U(0; x + e) = I. It is now con- 
venient to move to the interaction picture 



U(t;x + e) = U{t;x)W(t;x), 
which leads to a simpler equation 



(25) 



i — W(t;x) = eU(t:x^Hi(t;x)U(t;x)W(t;x). (26) 
at 

Integrating formally this differential equation and iterat- 
ing the resulting formula, one obtains the usual Dyson 
series. We only need this series up to first order 

W(t; x) = I - ie [ drU{T] x^H^t; x)U(t; x) + 0(e 2 ). 
Jo 

(27) 

From (f2"5")) combined with (f2"7) we have 



U{t;x + e) = U{t;x) - 



(28) 



ieU(t) / drf/(T) t ffi(r)J7(r)+0(e 2 ) 



(In the second term and hereafter we omit the x depen- 
dence to ease the notation). We can now compare this 
expression with (|22p in order to identify the derivative of 
the unitary, that is 

JLU(t) = -ieU(t) drU(T)^(T)U(r). (29) 

Therefore we obtain the formula for the gradient of the 
fidelity 



dF _ 1 
dxj d 



Re E [ dT (up{T)U{r)^{r)U{T, 



(30) 



which relates the gradient of F to a derivative of the uni- 
tary operator. What follows is a simple way to compute 
dU /dxi which is based on performing a Taylor expansion 
of the operator U(t; x) with respect to the parameters x 

BTI 

Ufa x + e) = U{t; x) + e—(t; x) + 0{e 2 ). (22) 
ox 

We will obtain this series using time-dependent pertur- 
bation theory on the Schrodinger equation (|18p . which 
will enable us to identify the derivative of the unitary 
operator. 

Let us assume that by changing x — > x+e our Hamilto- 
nian decomposes into an unperturbed part, Hq = H(t; x) 
and a perturbation H\ 

ri M 

effi = H(t; x + e)-H = e— (t; x) + 0{e 2 ). (23) 

The new unitary operator will satisfy a Schroodinger 
equation with a modified Hamiltonian 



C. Development of the algorithm 

Even though we have a closed expression for the deriva- 
tive of the fidelity, we still need to compute the integral 
which appears in Eq. (f30|) . We have devised a simple, 
accurate and efficient procedure which is based on solv- 
ing three sets of ordinary differential equations. The first 
one is obtained by transforming the integral in Eq. (|30p 
into d x M ordinary differential equation 



d_ 

df 



^Re (ip n 



Up(T)U(ty—(t)U(t) 



Vv 



with initial condition / n ,i(0) = and final value 



(31) 
(32) 



We now notice that left side of the scalar product in 
Eq. (EU is the state U(t)U(T)W g \ij; n ) and can be com- 
puted by solving an ordinary differential equation 



ijU{t- x + e) = (H + e#i)t/(t; x + e), (24) 



(33) 
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with initial condition \£ n (T)) — Ug\ip n ) and moving back- 
wards in time from T to t. 

We now have all the ingredients to design the protocol 
that computes our derivative (|30p . We summarize it as 
follows. First, solve Eq. (j33j) backwards in time up to 
t = 0, finding 

|£„(0)> := U{T)^U g \i> n ). (34) 

We then solve the system of ordinary differential equa- 
tions 

ij t \Mt)) ■= H(t)\Mt)) (35a) 

ij t \Ut)) ■= H(t)\Ut)) (35b) 

r] 1 f) M 

-/„,, := -I m (Ut)\^(t)\Mt))- (35c) 

With the initial conditions \ip n (0)) being some orthonor- 
mal basis of our Hubert space, |£n(0)) computed before 
and f n ,i(0) = 0. The value of the derivative is then com- 
puted using Eq. (|3"2")l . 

This derivative and the formulas given above can be fed 
to any optimization package, such as a simple line search 
algorithm or a nonlinear conjugate gradient method. In 
particular we have used Matlab's nonlinear optimization 
toolbox 

VII. CONCLUSIONS 

In this work we have proposed to generate quantum 
transport using cold atoms in an optical superlattice that 
is modulated periodically both in time and in space. 
Since there is no average force acting on the atoms, we 



deal with a quantum ratchet effect that makes atoms on 
even and odd sites move along opposite directions. 

We have demonstrated that the ratchet effect can be 
induced both for free particles and in the case of nonzero 
on-site interactions. The latter represents to our knowl- 
edge the first time strong many-body interactions have 
been treated exactly in Hamiltonian ratchets (5?| 

This ratchet effect can be used to transport quantum 
information by imprinting a qubit or an entangled state 
on one or more atoms of the atomic chain and letting the 
system evolve according to our ratchet potentials. The 
dynamics generated with our scheme corresponds to a 
smooth concatenation of perfect swap gates. Therefore, 
after a time n x T, we will find that the qubit state has 
been transported a well determined distance, of order n, 
along the lattice. 

Our ideas could be implemented in current experi- 
ments with optical superlattices. From the quantum in- 
formation point of view, compared to other proposals 
based on effective spin interactions between atoms, the 
ratchet effect should be faster and lead to a more flexi- 
ble dynamics. From the fundamental point of view, we 
believe that the modulated superlattices are a rich play- 
ground in which to study transport phenomena, quantum 
diffusion and the influence of noise and of chaos in the 
transport of quantum states. 
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